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Bursts of particle production during inflation provide a well-motivated mechanism for creating 
bump-like features in the primordial power spectrum. Current data constrains these features to be 
less than about 5% the size of the featureless primordial power spectrum at wavenumbers of about 
0.1 /i Mpc -1 . We forecast that the Planck cosmic microwave background experiment will be able 
to strengthen this constraint to the 0.5% level. We also predict that adding data from a square 
kilometer array (SKA) galaxy redshift survey would improve the constraint to about the O.f % level. 
For features at larger wave-numbers, Planck will be limited by Silk damping and foregrounds. While, 
SKA will be limited by non-linear effects. We forecast for a Cosmic Inflation Probe (CIP) galaxy 
redshift survey, similar constraints can be achieved up to about a wavenumber of 1.0 h Mpc -1 . 
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I. INTRODUCTION 



Current data is remarkably consistent with predictions 
of the simplest models of inflation. To a high degree 
of accuracy, the Universe appears to be flat and have 
nearly scale-invariant, isotropically distributed, Gaus- 
sian, and adiabatic primordial perturbations [TJ [5] . How- 
ever, there are still potentially large improvements to be 
implemented in the precision and the length scales of the 
primordial perturbations that will be probed. Therefore, 
it is important to investigate possible deviations from 
this simple picture that may in future be detectable. In 
this article, we concentrate on features in the primordial 
power spectrum that may be caused by bursts of particle 
production during inflation [3 5 . This can happen when 
the motion of the inflaton causes the mass of another field 
to pass through zero. The resultant production of par- 
ticles leads to a corresponding bump like feature in the 
primordial power spectrum at around the scales which 
are then leaving the Hubble horizon. 

A feature in the primordial power will translate to 
a corresponding feature in the cosmic microwave back- 
ground (CMB) angular power spectrum and the matter 
power spectrum. The matter power spectrum may be 
probed in many ways, and in this article we concentrate 
on galaxy redshift surveys and cluster surveys. 

Currently, there is no detection of such features in 
the data, but only wave-numbers of k < 0.1 h Mpc -1 
have been probed and only to accuracies of about 5% 
[3]. The Planck CMB survey can probe smaller length 
scales due to higher resolution and lower noise and so will 
help improve the constraints up to to k ~ 0.2 h Mpc -1 
where noise, beam size limitations and foregrounds start 
to dominate. In general astro-physical foregrounds pre- 
vent one probing the primordial power spectrum beyond 
about k ~ 0.5 /i Mpc -1 using the CMB. Future galaxy 
redshift surveys such as are planned with the Square Kilo- 



meter Array 1 (SKA), will also probe k ~ 0.1 ZiMpc -1 to 
much higher accuracy due to the huge volumes that they 
will encompass but will be limited to k < 0.2 h Mpc -1 by 
non-linear effects. While, k < 2 ft, Mpc -1 may be probed 
by very high redshift galaxy surveys, such as the Cosmic 
Inflation Probe 2 (CIP), where non- linear growth has yet 
to dominate. 

In this article, we will make forecasts on how well the 
amplitude of a particle production-induced feature can be 
constrained by future large scale structure surveys. We 
begin in Section [H] with a summary of how particle pro- 
duction can generate a bump-like feature. In Section [lll[ 
we give an overview of Fisher information matrices, a 
conventional method of predicting constraints on a set of 
parameters, for CMB, galaxy and cluster surveys. Pre- 
dictions for constraints on the feature amplitude and po- 
sition combined with other cosmological parameters are 
given in Section IV We discuss our results in Section [V] 



II. PARTICLE PRODUCTION 

Recently, a mechanism has been proposed that will 
generate a bump-like feature through particle production 
during inflation (3H2]. In this scenario, the production of 
massive iso-inflaton particles during inflation gives rise 
to potentially quantitatively observable features in the 
primordial power spectrum. The fields simply interact 
via the coupling 



(i) 



where g m is the interaction coupling constant. <p and \ 
are the inflaton and iso-inflaton fields respectively. When 
cj> passes through <j> there is a non-adiabatic change in 
the mass of x an d as a result a burst of particle produc- 
tion. This rapidly drains energy from the (f> field which 
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cluster number counts, the cluster/galaxy power spec- 
trum and the cosmic microwave background power spec- 
trum. 




1. Differential Halo Mass Function 
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FIG. 1: The primordial power spectrum plus particle 
creation features with Am = 1.25 xlO -10 at position 
fcm = 0.1/iMpc" 1 (red dashed). 



can lead to a transient violation of slow roll and hence an 
associated "ringing" in the primordial curvature fluctua- 
tions which is similar to that seen in models with a sharp 
feature in the potential, [6Tfl4]. However, the dominant 
effect is found to come from multiple re-scatterings of 
the produced S\ particles off the <j> condensate. Mul- 
tiple bump-like features scenario is also possible with 
</> 's in different positions, which associates with different 
<7 IR 's. However, we restrict to the case of only one feature 
presents. The overall effect is a bump-like feature which 
can approximately be described by a parametric form as 
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where is the scalar amplitude describing the normal- 
isation of the power spectrum. n s and fc p i vot is the tilt 
and the pivot wavenumber respectively. The first term 
on the RHS of Eq. ([2| is the standard power-law power 
spectrum and the second term contains the features gen- 
erated by particle creation which is parameterised by an 
amplitude, A m and position fc IR as shown in FIG.[T] The 
normalisation is chosen so that A IR is the amplitude of 
the feature at its peak fc IR . The relation to the coupling 
constant is given by 



A IU M 1.01 X 10 
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(3) 



III. FISHER MATRIX CALCULATION 

We use the Fisher Information Matrix [15] to make pre- 
dictions for constraints of cosmological parameters for fu- 
ture surveys. The statistics that will be implemented are 



In order to predict the number density of collapsed 
objects in the Universe, a statistical concept of halo or 
mass distribution is used here. The differential halo mass 
function, or halo mass function for short, is defined as 
the redshift-dependent distribution of the number of col- 
lapsed dark matter haloes per unit mass interval in a unit 
co-moving volume. The halo mass function is given by 



dn p m d In a 
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where dn/dM is the differential halo mass function, p m is 
the matter density, and /(c) is called the mass fraction. 
The smoothed variance is calculated as 



a 2 (R,z) = 



2tt 2 



P(k) W 2 (k,R) k 2 dk, (5) 



where D(z) is the linear growth function normalised to 1 
at the present epoch. W(k, R) is the Fourier-space top- 
hat window function. R is the smoothing radius for a 
comoving sphere enclosing a mass of 



M = %R*p„ 
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(6) 



The top hat smoothing in Eq. |5]) suppresses the con- 
tribution of any change to the primordial power spec- 
trum located at wave number fc IR 1/R. Combined 
with Eq. (((H), this implies that a change in the primor- 
dial power spectrum at fc IR has a suppressed effect on the 
number density on mass scales satisfying 
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(7) 



The mass fraction is defined as a fraction of mass in 
collapsed haloes per unit interval in lner . The halo 
mass function is described by a pair of parameters, f(o~) 
and lner -1 . Both of these parameters are a natural way 
of parameterising the mass function from different cos- 
mological models with the fewest number of parameters. 
All the cosmological parameters are embedded in a. f(a) 
is the fraction of the density fluctuation that eventually 
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collapse into non-linear objects. We used the mass func- 
tion by Jenkins et al. (2001) [16]: 



f(a) = 0.315exp [-| lncr" 1 + 0.61| 3 ' 8 ] . 



2. Number Count Fisher Matrix 



(8) 



For a survey which covers /sky fraction of the sky, a 
theoretically expected value of the number of clusters in 
the ith redshift bin at central redshift Zi and a width of 
Az is given by 



N(M>M lim ,z) = 

dV 

/sky / dz 



dz 



Ah 



lim 



dM^<M,;). 



The differential co-moving volume, dV^/dz, is 
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where -ff(z) is the Hubble parameter at redshift z. 
Eq. ([9| includes all the clusters above a threshold mass 
of Mii, n . The mass threshold can, in general, be redshift - 
dependent which is normally given in terms of a survey- 
specific selection function. To investigate generally how 
cluster surveys can be used to probe features, we will 
use a simple redshift-indcpcndcnt effective mass thresh- 
old for cluster surveys which yields an equivalent number 
of clusters over the entire survey volume. 

Given a set of parameters of interest, = 
(6 1 ,9 2 , ...,6 m ), an element M and 6„ of a Fisher matrix 
for number count is [17] 
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(11) 



where F^ u is a sum of all redshift bins in the survey. 



B. Power Spectrum 

We also consider cosmological constraints from mea- 
surements of the matter power spectrum. 



1. Galaxy Power Spectrum 

Galaxies are moving away from us along with the Hub- 
ble flow. The perturbation of the velocity of the galaxies 
is called the peculiar velocity which is the motion of the 
galaxies relative to the Hubble flow. It changes the ob- 
served Doppler shift and, hence, the distance inferred 
from redshift measurement. Kaiser (1987) [TB] showed 




(9) 



FIG. 2: The matter power spectrum from a featureless pri- 
mordial power spectrum (blue solid) and the matter power 
spectrum from a primordial power spectrum which has a par- 
ticle creation feature with amplitude Ai R = 1.25 xlO -10 at 
(10) position km = 0.1 h Mpc -1 (red dashed). 



that the power spectrum derived from a redshift survey, 
P s (k) , is given by 



P.(fc,M)= [l + f3 f i 2 ] 2 b 2 P(k), 



(12) 



where P(k) and P s (k) are the matter power spectrum 
and the redshift power spectrum respectively, fi = k • n 
is the cosine of the angle between k and the line of sight, 
P is 



_ ldlnD 
b din a 



(13) 



where D is the linear growth function normalised to 1 
at the present epoch, and b = 8 g /8 is the galaxy bias. 
The bias for galaxy surveys is normally estimated from 
the data and the uncertainty in bias is propagated into 
the other parameter constraints. However, to include 
non-linear effects, we applied a Taylor-like expansion to 
the bias as 



b 2 = bl [1 + a x k + a 2 k 2 ] 



(14) 



where bo is a scale independent bias. a 1 and a 2 are the 
first and the second order term respectively |19j . 

The linear theory matter power spectrum is related to 
the primordial power spectrum by 



P(k) oc T 2 (k) k V{k) 



(15) 



where T is the transfer function. So a narrow feature in 
the primordial power spectrum will be transferred to a 
narrow feature in the matter power spectrum at about 
the same comoving wavenumber, see FIG. [2] 
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2. Cluster Power Spectrum 



The Poisson approximation for cluster number counts 
is not strictly accurate for larger surveys 20] . There is 
an additional sample variance which can be used to help 
constrain the mass scaling relation. We account for this 
by incorporating an additional constraint from the power 
spectrum of the clusters. The way of doing this is similar 
to that of galaxies. However, in the case of cluster bias, 
we use an effective linear halo bias 



b e s(z) 



r~ dM b(M, z)dn/dM 

m lim 



J^dMdn/dM 



(16) 



where dn/dM is the differential mass function (See 
Eq. Q) and b(M, z) can be calculate from halo bias [5T] 



b{y) 



where v = (S c /a) 2 with 6 C 



1 



2p 



(17) 



5 c [l + (au)P}' 
1.69, a = 0.75 and p = 0.3 



3. Power Spectrum Fisher Matrix 

We may write the appropriate Fisher Matrix for the 
power spectrum by assuming the likelihood function to 
be Gaussian as [2H 



d 3 k <91nP s (k) TA ,ainP,(k) 
Voff (k) - 



(27T) 



(18) 



We set fc max to be the wavenumber where non-linear ef- 
fects start to become non- negligible. From 23- where 

We set 
The ef- 



c(i?ni) = 0.5 and k 
k min = 1.0 x 10~ 4 Mpc 
fective survey volume, V c g, is a given by 



!im:{ — k n \ — ir/2/R n \. 
1 for the lower limit. 



V cS (k,fi) = / d^r 



n(r)P s (k, /i) 



1 2 



n(r)P s (fc,/i) + 1 



nP s (fc, (i) 
nP,{k,fj) + 1 



Vs. 



(19) 



where n(r) is the number density at position r and n 
is the average number density for the survey. The ef- 
fective volume is due to a finite survey volume and in- 
complete sampling of the underlying density field. These 
are known as sample variance and shot noise respec- 
tively. The uncertainty is propagated through the cal- 
culation by the weighing factor [n(r)P s /(n(r)P s + l)] 2 . 
For galaxy surveys where the number density is high 
ie. nP s (k,(i) 3> 1, the effective survey volume is then 



eff 



V. 
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C. CMB 

The primordial power spectrum is probed over a wide 
range of wave numbers by measurements of the primary 




FIG. 3: The CMB angular power spectrum from a featureless 
primordial power spectrum (blue solid) and the CMB angu- 
lar power spectrum from a primordial power spectrum with 
a particle creation feature which has A IR — 1.25 xlCP 10 at 
position feia = 0.1/iMpc^ 1 (red dashed). 



CMB anisotropics (see for example [21] )• Both the tem- 
perature (T) and -E-mode of the polarization (E) probe 
scalar perturbations. 



£(£ + 1)C? 
2tt 



= J d\nk{T e x (k)) 2 V(k), 



(20) 



where X <ET,E for auto-correlation function (TT, EE) . 
For the cross-correlation power spectrum (TE), 



1(1 + l)Cf 
2tt 



d]nkT[(k)Tf(k)V(k). 



(21) 



The projection of a mode of wave-number k on to the 
surface of last scattering (a sphere of comoving radius 
D*) results in the CMB transfer functions having the 
form Tf - jt(kD*). Where ji is the spherical Bessel 
function of order t which peaks at £ ps kD* . Therefore a 
bump in the primordial power spectrum at wave-number 
fcj is mapped onto a bump in CMB angular power spec- 
trum at 



« 10 4 - 



(22) 



hMpc~ L 
see FIG. 1 

The foreground contribution from secondary sources 
which will probably be hard to completely remove for 
£ > 2000 for both temperature and polarization. For this 
reason, as done by [21] and [5S] , we will restrict ourselves 
to I < 2000 when evaluating the forecasted marginalized 
errors. 

The CMB Fisher matrix is given by (see for example 



.) 



F„ = 



i x,x< 



Gov" 1 (Cf,Cf ) 
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89, 



(23) 
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TABLE I: Planck Instrument Characteristics 



Center Frequency (GHz) 


70 


100 


143 


217 


9 (FWHM arcmin) 


14 


10 


7.1 


5.0 


<T T (/J.K) 


12.8 


6.8 


6.0 


13.1 


a E (a*K) 


18.2 


10.9 


11.4 


26.7 



where the covariance matrix can be obtained from [26] 
and it depends on the temperature noise per pixel (a T ), 
the polarization noise per pixel (er E ), the pixel area in 
radians squared (0 2 = 47r/A r p i x ), and the beam win- 
dow function which we approximate as Gaussian (B t ss 
exp(-^+lK 2 ). 



IV. FORECAST CONSTRAINTS 

With the descriptions of Fisher information matrices in 
Section [ill} we can make predictions for future upcoming 
surveys. We also test our Fisher matrix formalism on 
some current surveys and compare to published results. 
For the CMB we take the Planck 3 and WMAP 4 surveys. 
The values for Planck are taken from the Planck blue 
book 5 and are listed in TABLE [i] (note that 6 needs to 
be converted to radians). We use a b = 0/^8 log 2 and 
combine the different frequency bands as specified in |27j . 
We take the range in £ to be 2 to 2000. At higher £, 
secondary sources of temperature and polarization will 
likely prohibit the extraction of cosmological information 
from the primary CMB. 

For galaxy surveys, we consider the Cosmic Infla- 
tion Probe 6 (CIP) [25], which is a space-based mission 
aimed to measure the linear galaxy power spectrum over 
k ~ 0.03 - 2.0ft,Mpc -1 to better than 1%. The primary 
science goal of the CIP is to provide constraint on in- 
flation models by observing Ha. The CIP survey will 
cover 1,000 square degrees and will be capable of detect- 
ing more than 10 8 galaxies between redshift range of 1.8 
-6.5. We follow the CIP model from by having redshift 
bins at z — 2.0 — 3.5, 3.5 — 5.0, 5.0 — 6.5 and and hav- 
ing an average galaxy number density of 1.0 xl0~ 2 , 5.3 
xlO -3 and 1.3 xl0~ 3 h 3 Mpc -3 respectively which is 
equivalent to ^100 million of galaxies within the survey 
volume of ^15.0 ft,~ 3 Gpc 3 . We restrict our calculation 
to fc max = 0.5, 1.0, and 2.0/iMpc^ 1 for the three bins 
respectively (See Eq. (18l). 

The Sloan Digital Sky Survey 7 (SDSS) [2J is a ground- 
based optical survey using a 2.5 meter telescope at 
Apache Point Observatory, New Mexico. It final data 



3 http:/ /www.rssd.esa.int/index.php?projcct=planck 

4 http://map.gsfc.nasa.gov/ 

5 http:/ /www.rssd.esa.int/SA/Planck/docs/Blucbook-ESA- 
SCI(2005)l_V2.pdf 

6 See http://www.cfa.harvard.edu/cip/ 

7 See http: / /www. sdss.org/ 



set (DR7) includes more than 230 million celestial ob- 
jects and spectra of about 930,000 galaxies and 120,000 
quasars over an area of 8,400 square degree in five op- 
tical bandpasses to redshift about 1.0. We follow the 
SDSS survey model from Pritchard & Pierpaoli (2008) 
|30j . We estimate the bias for the SDSS survey as 2.25 
(ie. b = crg g/crg where a S g — 1.8 and er 8 = 0.8) and 
n = 1.0 x 10~ 4 ft, 3 Mpc" 3 . The Square Kilometre Ar- 
ray 8 (SKA) is a large-scale radio telescope which aims to 
cover the frequency range from 60 MHz to 35 GHz and 
20,000 square degree of the sky. The main science goal 
of the SKA project is to study the HI content of galaxies 
to cosmologically significant distances, z ~ 2.0, and to 
make reionisation maps using 21 cm transition of neu- 
tral hydrogen. The SKA project is currently in a design 
phrase and the telescope site will be decided in 2012. We 
consider the galaxy redshift survey component of SKA 
(see for example [31]). It also may be possible to use the 
21 cm absorption component of the SKA to probe reion- 
ization and matter power spectrum at a higher redshift 
(see for example [32]). We shall consider the SKA survey 
as a stereotype of a cosmic variance-limited galaxy red- 
shift survey by assuming that the number density is so 
high that nP s 3> 1.0. Hence, the effective survey volume 
is equivalent to the survey volume (see Eq. (19)) from 
z = 0.0 — 2.0. The set of parameter of all the galaxy 
surveys will include bias as an additional parameter re- 
flecting the fact that we do not know accurately the value 
of the bias. The uncertainty in bias will propagate into 
the parameter constraints. 

For cluster surveys, we consider an all-sky Extended 
ROcntgen Survey with an Imaging Telescope Array 9 
(eROSITA) which is a high sensitivity all-sky X-ray sur- 
vey in the 0.2-12 keV energy band. The key science goal 
for eROSITA is to constrain the properties of dark en- 
ergy using high redshift clusters of galaxies. It will have a 
capability to measure the spatial correlation features and 
evolution of a sample of about 50,000 galaxy clusters over 
a redshift range of 0.0 - 2.0 and will be able to find col- 
lapsed objects with mass above 3.5 x 10 14 h~ 1 M Q . Our 
calculation of the expected cluster counts is in a good 
agreement with the all-sky eROSITA cluster count give 
in eROSITA documentation. We model the eROSITA 
survey by having Az = 0.2 and z = 0.0 — 2.0. We also 
assume a 10% prior on the mass threshold determina- 
tion for the eROSITA survey. The set of parameters of 
eROSITA will include the mass threshold as an addi- 
tional parameter reflecting the fact that we do not know 
accurately the mass of a cluster. The uncertainty in mass 
threshold will propagate into the parameter constraints. 

The SuperNova/Acceleration Probe 10 (SNAP) survey 
[331 134) is a space-based experiment aiming to study dark 



8 See http://www.skatelescope.org/ 

9 See http:/ /www. mpe.mpg.dc/heg/www/Projects/EROSITA/main. html 
10 



See http:/ /snap. lbl.gov 
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energy and dark matter. We forecast the SNAP cluster 
lensing survey constraints [35] (also see [36H39] ) . Here we 
take z — 0.0 - 1.4, Az = 0.2, / sky = 0.024, and a mass 
limit of 10 li h~ 1 M Q . This roughly matches the number 
of clusters found using the more accurate selection func- 
tion of 35] when we use their fiducial model cosmological 
parameters. Also, [35] found that the signal-to-noise ra- 
tio of a more realistic selection function was about the 
same as taking a mass limit of 10 14 /i _1 M Q . Our SNAP 
selection function is biased to slightly higher redshifts 
than that of [35], but we expect this not to alter our 
predicted constraints significantly. We do not include 
the mass threshold into our analysis for the SNAP weak 
lensing survey as the weak lensing technique can be used 
to determined the cluster mass accurately. 

A summary of all galaxy and cluster survey parame- 
ters is given in TABLE [TTJ We combined both statistical 
measurements from cluster number counts and the power 
spectrum in our predictions for the eROSITA and SNAP 
surveys. For the cluster surveys, we estimated the effec- 
tive bias using Eq. (16). 



We follow the work done by Barnaby & Huang (2009) 
[3J where they found a 2-cr constraint on the ampli- 
tude of the feature of about 2.5 x 10~ 10 on scale k ~ 
0.006 — 0.1 /iMpc" 1 . We conservatively define our fidu- 
cial value for the amplitude as A m = 1.25 x 10~ 10 and 
the feature at positions fc IR = 0.1, 0.4 and 1.0 /iMpc -1 . 
A survey fc max limits the range of feature positions it 
can probe. The SDSS and SKA survey can only probe 
the feature at k = O.lftMpc -1 while CIP is the deep- 
est survey that can probe all of the above feature posi- 
tions. The eROSITA and SNAP survey probe only the 
k = 0.1/iMpc -1 and k = 0.4/iMpc -1 features, whereas 
the latter is probed by the number count component. 

We define our fiducial set of parameters as the flat 
ACDM set of parameters with running of the spectral 
index, a 3 , plus the amplitude and the position of the fea- 
ture O = (w b , u c , fl A , n s , a„ A|, r) = (0.0227, 0.1107, 
0.738, 0.969, 0.0, 2.15 x 10~ 9 , 0.086) and k pivot = 0.05 
Mpc -1 . Our fiducial cosmological parameters are con- 
sistent with the WMAP7 maximum likelihood values [2] 
although their fc p i vo t = 0.002 Mpc -1 . We assume a neg- 
ligible neutrino mass, but we checked that adding a non- 
negligible mass would not significantly change our results. 
The constraints for the fiducial parameters without a fea- 
ture are given in TABLE |III| For features, we show the 
marginalized constraints on A 1R , k m , n, and a s in TA- 
BLE |lVj TABLE |Vj and TABLE [VI] for features at fc IR = 
0.1, 0.4 and 1.0 ft. Mpc -1 respectively. The other cosmo- 
logical parameters are included in the marginalization, 
but they are not particularly degenerate with the fea- 
tures. We use a modified version of the CAMB package 11 
[40) to generate the power spectrum for our analysis. 
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FIG. 4: Marginalised probability contours containing 68% of 
the posterior probability between ot s and Am = 1.25 x 10 -10 . 
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k m [h Mpc ] 

FIG. 5: l-o" marginalised constraints on the amplitude and 
position of particle production feature for Planck, Planck + 
SKA and Planck + CIP. The uncertainty in fci R is multiplied 
by 5 to make them more visible. 



V. DISCUSSION 

For the standard cosmology, our cosmological con- 
straints are in good agreement with previous work. For 
example, our CIP constraints for n s and a s are consistent 
with [35] [3T], our results for SDSS and SKA are consis- 
tent with 30 . Our best constraint for (n„,a s ) derives 
from CIP + Planck and SKA + Planck which are (0.0020, 
0.0032) and (0.0016, 0.0019) respectively. However, since 
CIP and SKA surveys probe the power spectrum at ex- 
clusively different redshift range, we can consider them 
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as two independent surveys. The combined CIP + SKA 
+ Planck improves the constraints to (0.0014, 0.0017). 



From TABLE [TV] - |VIJ we investigate the effect of a 
feature on cosmological parameters on scales k = 0.1, 0.4 
and 1.0 /iMpc -1 respectively. The WMAP+SDSS con- 
straints for the O.l/iMpc -1 feature are consistent with 
those found for the actual data [3]. As the CMB is 
limited to I < 2000, it can only directly constrain the 
k = O.l/iMpc -1 feature. It may be possible to ex- 
tend this range somewhat for polarization, but proba- 
bly not enough to completely encompass the width of 
the 0.4 /iMpc -1 centred feature. The cluster surveys 
were only able to help directly constrain the feature at 
0.1 /iMpc" 1 . Although in previous work [42] . we found 
that cluster surveys could constrain the primordial power 
spectrum directly at scales of 0.4/iMpc -1 , this was for 
a linear piecewise binning of the primordial power spec- 
trum where the location of the feature did not need to 
be constrained. 

With features at different positions, the constraint on 
n, remains almost the same while the constraints on a s 
are worsened by about a factor of 2. The degeneracies 
between A 1R and k ln have a significant impact on the 
running as a variation in a„ could be made on a more 
localised scale than a variation in n s . However, there 
is more of an improvement on a B with CIP + Planck. 
FIG. |4] shows that CIP could be used to break the de- 
generacy between A m and a s while SKA only provides 
a tightened constraint but the degeneracy is still in the 
same direction. CIP is in a different direction as it has 
most of its constraining power at higher k values and so 
the degeneracy between A m , A^, and a s is different. 

Our inference of a A ,„ and er,. TT , is summarised in FIG.[5l 

^IR ^IR I I 

As can be seen, Planck will improve the constraint at the 
0.1 /iMpc -1 scale by about a factor of about 10 in com- 
parison to WMAP. The addition of SKA further improves 
the constrain by another factor of about 5. On wavenum- 
bers up to about 1.0 h Mpc -1 CIP combined with Planck 
could provide similar constraints. 

Even smaller scales are constrainable, but then CIP 
would not be able to probe the full extent of the fea- 
ture. Also, it may be possible to further improve the 
constraints by including the non-Gaussianity associated 
with the feature [5]. 

Even though the excess in cluster number count can 
be easily obtained from cluster surveys (See FIG. [6] for 
eROSITA), they do not give better constraints in com- 
parison to galaxies surveys partly due to low number 
statistics. In addition, the cluster number counts are not 
good at simultaneous determining the amplitude of the 
features and the positions due to the degeneracy in clus- 
ter counts. Because the number of clusters is determined 
by the integral of P(k) (See Eq. Q and Eq. ^) the fea- 
tures with different amplitude and position can conspire 
to yield the same integral, hence, number counts. FIG. [7] 
shows the contour plot of excesses in number counts for 
different amplitude and position of the features. 

In summary, we have demonstrated that future surveys 



• • No Feature 
» t k m =0.1hMpc" 
■ ■ fc Ia =0.4/iMpc~ 




FIG. 6: The effect of a change in primordial power spectrum 
on the eROSITA all-sky survey. The circle dots represent 
number of cluster that would be observed by eROSITA for 
no feature power spectrum. The triangle and square dots 
show excess in number count for feature at km — 0.1 and 0.4 
ftMpc" 1 respectively. The error bars are la. 




FIG. 7: The excess number count from SNAP survey due to 
the feature with an ampitude Air and position at fc m . 



can potentially probe the primordial power spectrum, for 
particle production-induced features, significantly more 
accurately and to significantly smaller length scales than 
at present. Even if no features are detected, at least 
the simplest models will have been tested much more 
precisely and on a much greater range of scales. 
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TABLE II: Details for galaxy and cluster surveys. 



Survey 


z 


A^bin 


T/ 

v survey 


"'max 





M \im 


J sky 








(h~ 3 Gvc 3 ) 


( hMvc^ 1 ) 




(h~ 1 M Q ) 




Galaxy Survey 
















CIP 


2.0 - 6.5 


1.5 


15.0 


2.0 


1.0 


N/A 


0.024 


SDSS 


0.0 - 0.6 


0.6 


1.0 


0.1 


2.25 


N/A 


0.3 


SKA 


0.0 - 2.0 


2.0 


100 


0.4 


1.0 


N/A 


0.5 


Cluster Survey 










(6eff) 






SNAP 


0.0 - 1.4 


0.2 


2.76 


0.15 


2.0 


1.0 xlO 14 


0.024 


eROSITA 


0.0 - 2.0 


0.2 


220 


0.2 


9.0 


3.5 xlO 14 


1.0 



TABLE III: Constraints on fiducial cosmology 



Survey 










ffo s 


(xlO" 9 ) 


0V 


O M lim 

(xlO^r'Mo) 




(h- 1 Mpc) 


(h- 2 Mpc 2 ) 


Galaxy Survey 
























CIP 


0.0022 


0.0086 


0.0043 


0.027 


0.0082 


0.23 


N/A 


N/A 


0.0057 


0.039 


0.024 


CIP + WMAP 


0.00026 


0.00070 


0.0015 


0.0047 


0.0037 


0.024 


0.0068 


N/A 


0.0055 


0.037 


0.021 


CIP + Planck 


0.000099 


0.00026 


0.0012 


0.0020 


0.0032 


0.013 


0.0030 


N/A 


0.0043 


0.029 


0.019 


SDSS + WMAP 


0.00046 


0.0036 


0.020 


0.021 


0.017 


0.055 


0.012 


N/A 


0.17 


4.44 


32.89 


SDSS + Planck 


0.00013 


0.0010 


0.0055 


0.0031 


0.0049 


0.020 


0.0044 


N/A 


0.15 


4.33 


32.82 


SKA 


0.00064 


0.0028 


0.0016 


0.0096 


0.0033 


0.077 


N/A 


N/A 


0.0029 


0.028 


0.041 


SKA + WMAP 


0.00017 


0.00061 


0.00055 


0.0028 


0.0020 


0.018 


0.0055 


N/A 


0.0025 


0.026 


0.037 


SKA + Planck 


0.000082 


0.00016 


0.00035 


0.0016 


0.0019 


0.0069 


0.0018 


N/A 


0.0022 


0.022 


0.032 


Cluster Survey 
























eROSITA + WMAP 


0.00037 


0.0021 


0.011 


0.013 


0.0082 


0.046 


0.011 


0.19 


N/A 


N/A 


N/A 


eROSITA + Planck 


0.00013 


0.001 


0.0052 


0.003 


0.0041 


0.020 


0.0043 


0.09 


N/A 


N/A 


N/A 


SNAP + WMAP 


0.00039 


0.0017 


0.010 


0.017 


0.015 


0.042 


0.010 


N/A 


N/A 


N/A 


N/A 


SNAP + Planck 


0.00012 


0.00066 


0.0036 


0.0027 


0.0049 


0.018 


0.0042 


N/A 


N/A 


N/A 


N/A 


CMB 
























WMAP 


0.0006 


0.0057 


0.032 


0.030 


0.021 


0.06 


0.012 


N/A 


N/A 


N/A 


N/A 


Planck 


0.00013 


0.0011 


0.0057 


0.0031 


0.0049 


0.02 


0.004 


N/A 


N/A 


N/A 


N/A 



TABLE IV: Forecasted 1-cr marginalised uncertainties for 
fiducial k m = 0.1 h Mpc -1 , A m = 0.125 x 1(T 9 . 



Survey 


^IR 




0~n s 


(7 Qs 




(xlO^ 9 ) 


(ftMpc^ 1 ) 






Grcilcixy Survey 












n m 


U.UUo ( 


n no7 


U.UUoo 


pip i \X7"\/r A D 


n nno c: 
U.UUoO 


n nno a 


U.UU4 f 


U.UUoo 


v^. ±1 -\- i IcLLILJy 


u.uuuu 


o nni ^ 


0094 




srms -U WMAP 


n ifi 
u. io 


n 097 


Ofi^ 

U.UUO 


U.UO'i 


SD^S -U PlanrV 
OL'uJ 1 IdllCK 


01 3 

U.UIO 


n nnoQ 
u.uuzo 


OOzH 


n nn7Q 

u.uu i y 


9"K" A 

O I \ - V 


u.uuoy 


U.UU10 


n nns^ 

U.UUoo 


u.uuoy 


QKA 4- WMAP 


n 00^9 


n nni A 


00*30 
u.uuou 


n nn9R 

U.UUiU 


SKA -(- Planrk 


n ooso 

U . UUOU 


0.0010 


0.0016 


n nn93 


V^lUbLcI O ILL Vcy 










r-RO^TTA _i_ WA/TAP 


n in 

U. 1U 


n 09^ 
u.uzo 


0*33 

u.uoo 


01^ 
U.UIO 


eROSITA + Planck 


0.011 


0.0022 


0.0040 


0.0057 


SNAP + WMAP 


0.075 


0.027 


0.038 


0.025 


SNAP + Planck 


0.012 


0.0023 


0.0038 


0.0075 


CMB 










WMAP 


0.19 


0.027 


0.064 


0.035 


Planck 


0.013 


0.0023 


0.0044 


0.0079 



TABLE V: Forecasted l-cr marginalised uncertainties for fidu- 
cial fci R = 0.4 h Mpc -1 , A m = 0.125 x 10~ 9 . 



Survey 


(xlO" 9 ) 


(ftMpc^ 1 ) 






Galaxy Survey 
CIP 

CIP + WMAP 
CIP + Planck 


0.017 
0.0075 
0.0072 


0.015 
0.0070 
0.0066 


0.026 
0.0048 
0.0021 


0.0084 
0.0039 
0.0033 



TABLE VI: Forecasted 1-cr marginalised uncertainties for 
fiducial k m = 1.0 h Mpc -1 , Air = 0.125 x 10~ 9 . 



Survey 












(xlO" 9 ) 


(ftMpc^ 1 ) 






Galaxy Survey 










CIP 


0.031 


0.096 


0.027 


0.0096 


CIP + WMAP 


0.023 


0.087 


0.0046 


0.0043 


CIP + Planck 


0.022 


0.086 


0.0020 


0.0036 



